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Abstract 

We study linear and nonlinear wave propagation in a dense and cold hadron gas and also in 
a cold quark gluon plasma, taking viscosity into account and using the Navier-Stokes equation. 
The equation of state of the hadronic phase is derived from the nonlinear Walecka model in the 
mean field approximation. The quark gluon plasma phase is described by the MIT equation of 
state. We show that in a hadron gas viscosity strongly damps wave propagation and also hinders 
shock wave formation. This marked difference between the two phases may have phenomenological 
consequences and lead to new QGP signatures. 
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I. INTRODUCTION 



Before the appearance of the RHIC data, a tremendous effort had been dedicated to the 
search of quark gluon plasma (QGP) signatures and to determine ways of experimentally 
disentangling the QGP from a hadron gas (HG) Among the most promising signatures 
we could find J / ip suppression and enhancement of strangeness production. All the proposed 
signatures were based on medium properties which were very different in each one of the 
two phases. Charmonium suppression, for example, was based on the QGP properties of 
deconfinement and plasma screening [2J . Strangeness enhancement was based on the chiral 
symmetry restoration, which is fully realized in QGP but only partially in a hadron gas 

In the most recent years the search for QGP signatures had lower priority, mostly because 
we believe that QGP has been already produced at RHIC and at LHC and now the focus 
should be on the study of its properties ^|. We believe that it is still interesting to look 
for QGP signatures and in this work we take the first steps to propose a new one, based on 
viscosity. 

Recently we have learned many things about the QGP and one of the most striking is that 
it is an almost perfect fluid with very small viscosity [4|, |5[ . We have also learned, albeit with 
much larger uncertainties, that the hadron gas has a very large viscosity y, |7| . According to 
some estimates, the viscosity coefficients (both bulk and shear) in the two phases can differ 
by two orders of magnitude or more! This difference is so big that it motivates us to look 
for observables which are sensitive to it. So far the existing calculations have only addressed 
the effects of viscosity in the QGP phase and the shear and bulk viscosity coefficients were 
varied in a limited range. These changes produce some visible but not very large effects on 
global observables such as rapidity and transverse momentum distributions, flow coefficients 
and particle correlations js]. 

We will investigate the viscosity effects on wave propagation in a quark gluon plasma 
and also in a hadron gas. These waves may be caused, for example, by fluctuations in 
baryon number or energy density. These fluctuations may be produced by inhomogeneous 
initial conditions, which, as pointed out in are the result of quantum fluctuations in the 



densities of the two colliding nuclei and also in the energy deposition mechanism. T 



ncsc 



fluctuations and their phenomenological implications have been studied extensively 
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because they may be responsible for the angular correlations of particle emission observed in 
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heavy- ion experiments. There are also hydrodynamic fluctuations |q|, which are the result of 
finite particle number effects in a given fluid cell. This generates local thermal fluctuations of 
the energy density (and flow velocity) which propagate throughout the fluid. Furthermore, 
there may be fluctuations induced by energetic partons, which have been scattered in the 
initial collision of the two nuclei and propagate through the medium, loosing energy and 
acting as a source term for the hydrodynamical equations. Finally, there may be also freeze- 
out fluctuations, which may be caused by finite particle number effects during and after the 
freeze-out of the hydrodynamically expanding fluid. 

According to our current picture of high energy heavy ion collisions, a dense and hot 
hadron gas may exist at later stages of the collisions, when the fluid has cooled down and 
the quark-hadron phase transition has taken place. A dense hadron gas may also exist and 
play a more important role in lower energy collisions, such as those which will be performed 
at FAIR [jjj] . Finally, a dense and cold hadron gas is expected to exist in the core of compact 
stars. It is not clear which mechanism might be responsible for the production of waves in 
the hadronic phase, but it is nevertheless necessary to know how waves propagate in a HG. 
In particular, we need to know to what extent waves generated in the QGP survive the 
passage through a hadron gas. 

At intermediate energy proton-nucleus reactions, where we do not expect to form QGP, 
the impinging proton may be absorbed by the nuclear medium creating a perturbation in 
baryon density, which then propagates through the target. This wave may have a relatively 
large amplitude and must then be treated beyond the linearization approximation. This was 
first done in 16], where it was found that, for a perfect fluid and for a particular equation of 
state, a Korteweg - de Vries (KdV) soliton is formed. In jl^l a similar analysis was carried 
out with an equation of state derived from the nonlinear Walecka model with the same 



conclusion. The formalism was subsequently extended to relativistic hydrodynamics in 18], 
to the case of spherical waves in 19| and finally to a hadron gas at finite temperature in 20], 
where a more detailed numerical solution of the differential equations was performed and a 
wider class of relativistic mean field models was investigated. For some equations of state 
and/or some approximation schemes we find KdV solitons and for others we do not find 
them, but even when there is no soliton we always observe that localized pulses pr opag ate 
for long distances (from 10 up to 20 fm !) without loosing strength. In 21| and 22| we 
completed the studies of perfect fluids including the quark gluon plasma at zero and finite 
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temperatures. Even for the QGP, a fluid with a very different equation of state, the main 
features of pulse propagation were still very similar to those found in the previous works. 



A really strong difference was observed in 23[ where viscosity was introduced for the first 



time in our study of waves. In that work we studied the time evolution of cylindrical "hot 
spots" (or simply "tubes") in a hot and viscous quark gluon plasma. We could clearly see 
that a small viscosity strongly damps and broadens the tubes during their expansion and 
they are more easily mixed with the background fluid. 

In the present work we perform a systematic comparison between waves (both linear 
and nonlinear) in a hadron gas and in a quark gluon plasma, focusing on the effects of 
viscosity. For simplicity we shall consider only the case where the temperature is zero and 
derive the wave equations for perturbations in a non-relativistic viscous fluid. We will then 
find the solutions of these equations and see how they change when we change the viscosity 
coefficients. 

In the framework of non-relativistic fluid dynamics we can derive a theory of dissipative 
fluids which is remarkably successful: the Navier-Stokes (NS) theory 24j-|26[. We can use 
this theory to investigate the evolution of density perturbations in a hadron gas and in 
a quark gluon plasma. Perturbations are usually studied with the linearization formalism 
0j which is the simplest way to study small deviations from equilibrium and to obtain 
a wave equations. The propagation of perturbations through a QGP has been investigated 
in several works with the help of a linearized version of the hydrodynamics of perfect fluids 
10l . Ill] and also of viscous fluids 12[|. We will later go beyond linearization and study the 
effects of viscosity on the propagation of nonlinear waves. To this end we employ the well 
established reductive perturbation method (RPM) |28l-l33|. 

This work is organized as follows. In the next section we present the non-relativistic 
viscous fluid dynamics. In section III we briefly review the main expressions of the equation 
of state of both the hadron gas and the quark gluon plasma at zero temperature. In section 
IV we combine the hydrodynamical equations with the equations of state and derive the wave 
equations and present a numerical study. The last section is devoted to some concluding 
remarks. 
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II. NON-RELATIVISTIC VISCOUS FLUID DYNAMICS 



In what follows we employ the natural units c = 1 and H — 1. The metric is g 



fW 



diag(+, — , — , — ). The non-relativist ic Navier-Stokes equation is given by 



2im\ 



dv 1 , dv 1 1 dp 1 dli k% , . 
1- V K — £_ (l) 

dt dx k pdx 1 p dx k 



with the viscous tensor given by: 



n^-,fS + ^-|^-C^ (2) 



\dx k dx 1 3 dx 1 J dx 1 
The viscous coefficients are rj for the shear viscosity and ( for the bulk viscosity. We shall 
consider (i — k — I — x) for the one dimensional cartesian case, so the Navier-Stokes 
equation (JTJ) with fl2]) becomes: 

dv x dv x 1 dp If 4 \ d 2 u 



dt x dx pdx p\ 3 / dx 2 
where p is the mass density The continuity equation for the mass density p is given by: 

^ + V-( P v) = (4) 

In the nuclear medium, the mass density is related to the baryon density through p = Mps, 
where M is the nucleon mass. The continuity equation for the baryon density in one cartesian 
dimension is: 

dpB v dpB_ dvx_ __ q 
dt x dx dx 

III. EQUATION OF STATE 
A. Hadron gas 

The theory of nuclear matter has experienced a continuous progress. Recent developments 



include the use of effective field theory with chiral power counting (see for example 34J). 
For our purposes it is enough to work with a theory which reproduces the main features 
of dense nuclear matter. In this section we review the derivation of the equation of state 
from the nonlinear Walecka model 



35 



37| of cold nuclear matter. The Lagrangian density 



is given by: 

C = ^(id* - g v Vn — (M — g s <p)}^ + X - (d^d^ - m 5 2 2 ) - \f, v F^+ 



(6) 



where F^ u = d^V u — d u V^. The degrees of freedom are the baryon field if), the neutral 
scalar meson field <f> and the neutral vector meson field V^, with the respective couplings 
and masses. The equation of state is obtained through the usual mean field theory (MFT) 
approximation which consists in considering the meson fields as classical fields. This classical 
approach is based on the assumption that the baryonic sources are intense, their coupling 
to the meson fields is strong and the infinite nuclear matter is static, homogeneous in space 
and isotropic. From the calculations performed in 
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2m v 2 



2 rn s 2 , %r ,« N2 AM-M*f (M-M, 
(M-M*) 2 + b- — - ; +c- ; — 



we have the energy density: 

kp 



r*\4 



PB + 



7s 



2g& 



3g& 



^9i 



(27T) 



d 6 k \ k 2 



M* 
(7) 



where 7 S = 4 is the nucleon degeneracy factor and the effective mass of the nucleon is defined 



as M* = M — gs 4>q. In numerical calculations we will use 35 



37j M = 939 MeV, g v = 9.182 



and my = 783 MeV . The baryon density (pb) up to the Fermi level (kp) is given by: 



ip^ip 



PB 



Is 3 



6tt 2 



(8) 



The nucleon effective mass is determined by the self-consistency relation obtained from the 
minimization of e(M*) with respect to M*\ 



M* 



M — 



9s Is 
m s 2 (2tt) 3 



d A k 



M* 



+ 



9s 



Vk 2 + M* 2 m s 



*\3 



9s^ 



(M - M*Y + —AM - M*) 
9s 4 



(9) 



B. QGP 



The equation of state of the plasma phase is derived from the MIT bag model which 
describes a gas of noninteracting quarks and gluons and takes into account nonperturbative 
effects through the bag constant B 32(. This constant is interpreted as the energy needed 
to create a bubble (or bag) in the QCD physical vacuum. The baryon density is given by: 



1 1Q 

PB = 



3 (2tt) : 



d 3 k [n 



k n k\ 



(10) 



with 



n k = n k( T ) 



1 

1 _|_ e {k-\ti/T 



and 



n k = H( T ) 



l 

1 + e (k+^)/T 



(11) 
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where /i is the baryon chemical potential. At zero temperature the expression for the baryon 
density reduces to: 

P B = g^A*" (12) 

and h,F refers to highest occupied momentum level. The energy density and the pressure are 
given respectively by: 

e = B + ^3 J d 3 k k {e k ' T - I)" 1 + ^y 3 jd 3 kk [n s + %] (13) 

and 

p = ~ B + H / d " kk ( efc/T - l ^ 1 + j^fl d ' kk b + ^] } (14) 

The statistical factors are 7^ = 2 (polarizations) x 8 (colors) = 16 for gluons and 7q = 
2(spins) x 2(flavors) x 3(colors) = 12 for quarks. From the above expressions we extract the 
speed of sound c s : 



For the cold plasma at zero temperature and high baryon density the quark distribution 
functions become step functions. Inserting (fT2"j) into (fT3"]) and into (fT4j) we find the energy 
density and pressure as functions of the baryon density, respectively: 

e( PB ) = g) 7/ >W /3 + B (16) 

and 

p(Pb) = 1(1) V \ 2/3 Pb 4/3 -B (17) 

IV. WAVE EQUATIONS 
A. Linear waves 

Before addressing the nonlinear wave equations, we revisit the linear wave equation ap- 
proach in viscous hydrodynamics. The non-relativistic Navier-Stokes equation may provide 
a simple wave equation for perturbations in the baryon density based on the linearization 



formalism 



24 



27] . The use of this formalism can be justified when small perturbations are 
considered in a fluid at rest with background pressure po and background baryon density 
Pq. For simplicity we consider perturbations which depend only on the x space coordinate 
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and on time. We first define the dimensionless variables for the baryon density, for the fluid 
velocity and for the pressure: 

p(x,t) = ^M , ^ (M ) = ^) and = (is) 

Po c s p 

The expression for the speed of sound, c s , is: 

c s 2 = ^ p~c 2 e (19) 

In the non-relativistic limit e = p, where p is the matter density, which is related to the 
baryon density through p = Mps- We can then write: 

p = c s 2 Mp B (20) 

Rewriting ([T]), flS} and ( 1211 in one cartesian dimension and using f[T8~j) we find: 

„, J dv x 2 „ dv x \ dp / 4 \ a 2 ^ 

MpoP^— + c s t , x _j+ P0 _=^ + - l7 j c ._ (21) 

<9/5 A (9/5 A 9?) z „ 

i +CA i +c -"i^ =0 (22) 

and 

PoP = Mp c s 2 p (23) 
The perturbations are described by the following expansions of the dimensionless variables: 

p = 1 + <5/Ob , v x = 5v x and p = 1 + 5p (24) 

where <5ps, 5f x and Sp, denote small deviation of the baryon density, velocity and pressure 
from their equilibrium values respectively. When the expansions 024p are inserted into 
the relevant equations in consideration, the "linearization approximation" is performed by 



23| . Inserting ([24]) into (J2U), 022), flU and performing the 



neglecting the 0{5 2 ) terms 
linearization, we find respectively: 

d d ( 4 \ d 2 

Mp c s —Sv x + p —5p= {C + ^V)Csg^5v x (25) 

d d 

q1 5 Pb + c 'fa Sv * = ( 26 ) 

and 

p (l + 5p) = Mp c 2 (l + 5p B ) (27) 



Differentiating (|26|) with respect to x and inserting it into (!25l) we find: 

Mp c,|^ + Po ^-5p = - (C + ^) ||Up B (28) 
Performing the derivative of ( 1281) with respect to x and using ( 1271) , we find: 

M "» 41^ + M "° c ' 2 J^* = - ( c + 5") S5?*«> (29) 



Inserting the time derivative of (j26|) into 



for the baryonic density perturbation 



f29]) we obtain the following viscous wave equation 



38 



6 PB - -^l^tiPB = -V^TT^SPB (30) 



dx 2 c 2 dt 2 dt dx 2 

where 

and making use of the Ansatz: 

5p B (x,t) =Ae ikx ~ iwt (32) 
in (130]) we find the following dispersion relation u(k): 

u 2 = c s 2 k 2 — ivuk 2 (33) 

where v > and fcGR. The angular frequency is decomposed in two components j38|: 

oj = ujr + iuji (34) 

where Ur G R and U/Gi. Inserting (l34"j) in ( )33|) we find: 

cur = c s k — and ujj = — < (35) 

From the above we draw two important conclusions concerning the velocities: 
i) The phase velocity v p is given by: 

uj r \ 2 2 ( vk\ 



We can observe that one of the effects of viscosity is to reduce v p . 
ii) The group velocity v g is calculated using \u\ = c s k : 

d\uj\ 



Va = ~dk = ° s (3?) 



9 



Inserting ( |35|) into (134]) and the resulting value of oj into the Ansatz (|32|) . we find the solution 
of dSQ]): 

5 PB (X, t) = Ae -^t/2 ei (k X -t^c s ^-u^/,) (3g) 

with 

ft{5p B (x,t)} = Ae- vkHl2 cos(kx - \Jc s 2 k 2 - ^-tj (39) 

Since in the derivation of the above equation we have only used the general expression (fi~9"]) . 
it is valid both for the QGP and HG phases. The differences appear because the constants v 
and c s are different. The former varies over a much wider range of values whereas the latter 
is never very different from 1/3. For the sake of simplicity we will choose a constant speed 
of sound c s 2 = 1/3 for the two phases and all the differences will come from v. From (13 ip we 
see that the changes in density and viscosity are correlated in such a way that in a hadron 
gas r\ and ( are large and po is relatively small with v being large, whereas in the plasma 
phase the opposite happens and v is small. As an example we study the effects of viscosity 
on the waves propagating in a cold and dense system. In this case the perturbation 5p B is 
a deviation of the baryon density from the background density po- 

In numerical calculations the viscosity coefficients £ and r] may be taken from js, 7] and 
used in (13TT) . Using the exact values of these coefficients is not really crucial for our purposes. 
It is enough to know that rj is dramatically different in different phases. The phase velocity 
must be real and therefore in (139]) we must have: 

< — (40) 

From this inequality we see that, since vqgp << ^hg the maximum value of k is much 
smaller in the HG than in the QGP and consequently in a hadron gas there are much less 
modes than in a quark gluon plasma. From ( 139]) we see that for a hadron gas with high 
viscosity there is a strong damping of the wave due to the exponential factor, whereas in 
the plasma the oscillations persist for a very long time. 

We close this section emphasizing that waves produced by small perturbations, those 
which can be treated with the linearization approximation, behave quite differently in a 
hadron gas and in a quark gluon plasma. This pronounced difference comes essentially from 
the very different viscosity coefficients in these two phases. In a hadron gas waves are damped 
and stalled, whereas in a quark gluon plasma they can propagate nearly undisturbed. 
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Since there are many possible sources of perturbations and since they do not have to 
be necessarily very small, this study has to be extended to stronger perturbations and the 
resulting nonlinear waves. This is the subject of the next section. 



B. Nonlinear waves in a hadron gas 

We start recalling the relation between the mass density and baryon density, p = Mps, 
which used in equation (j3]) yields: 

( dv x dv x \ 1 dp 1 /. 4 \ d 2 v x 
Pb {^ + V ^)=-M^ + m[ C+ 3 71 )^ (41) 

From the first law of thermodynamics at zero temperature we have: 

de 

de = UbcIob where ug = tz — (42) 

OpB 

Inserting ( |42p into the Gibbs relation at zero temperature: 

de + dp — p B dp B + Pb dps 

we find 



dp = p B dp B 



and then finally: 



dp = pBd^^j and consequently ^ = p B ^ {^p~~^j 

Inserting f )43p into (141 j) we find: 

(dv x dv x \ 1 d ( de \ 1 / ^4 \d 2 v x 

where the relevant quantity which carries information about the equation of state is de/dps, 
which can be computed from (J7|). In the present calculation we will neglect the weak 
dependence of M* and also of the integral term in ((7|) on ps- We obtain: 

de g v 2 



dp B m v 2 

Inserting (l45p into (jUJ) we find: 



Pb (45) 



/ dv x dv x \ 1 



gv 2 dp B f A\ d 2 v 



my 2 dx \ 3 / dx 2 
11 



(46) 



which is the Navier-Stokes equation for the hadron phase. 

We next combine (1461) and (jSJ) to obtain a nonlinear wave e quat ion. In dealing with these 
equations we use the reductive perturbation method (RPM) [28M33I]. which is a technique 



which allows us to go beyond linearization and to preserve nonlinearities, dissipative and 
dispersive terms in the wave equations. We have already used the RPM in our previous works 
to study nonlinear waves in relativistic and non-relativistic ideal hydrodynamics and also 



to study the evolution of flux in relativistic viscous hydrodynamics 



23] . A comprehensive 



review of the technique and its applications in hydrodynamics of strongly interacting fluids 



can be found in 



32| 



The background baryon density is p , upon which the perturbation propagates. We now 
rewrite fl5J) and (JIB")) using ([IS]) to obtain: 



dp „ dp n dv x 
-rrr + c a v x — + c s p— 
at ax ox 







and 



Following 



P 



28 



dv x 
at 



+ c s v x 



dvx 
dx 



Po ~( 9v 
P 



Mc 



dp 
2 I dx 



3l|, we define the "stretched coordinates' 1 
£ = a V ^-c s t) and 



Mp 



T 



d 2 v x 
dx 2 



L 



(47) 



(48) 



(49) 



where L is a characteristic length scale of the problem (typically the radius of a heavy 
ion) and a is a small (0 < a < 1) expansion parameter. We also perform the following 



transformation of the viscosity coefficients 



39 



40) 



c 



a 



1/2 



and also 



(50) 



From the stretched coordinates (149]) we have the operators: 
a 1 ' 2 d d 2 a d 2 



d_ 

dx 



and 



d_ 
at 



a l/2 q a 3/2 q 

c x -^— + c s — r - — (51) 



L d£ ' dx 2 L 2 d£ 2 at L <9£ L dr 

Using (150]) and (]5T]) we transport f)4T]) and (HE]) from the (x, t) space to the (£, r) space, 
obtaining: 



dp 
d£ 



dp . . dp „dv x 
a a^ + V ^ + P ^ = ° 



and 



dv x 
d£ 



a 



dv x 
dr 



+ v x - 



dV-y 



d£ 



9v 2 po 



dp 
PlB + 



a 



Mc 2 m v 2 ) r di ' Mp c s L\^ + 3^ J d£ 2 



d 2 v. 



(52) 



(53) 
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We now perform the expansion of the dimensionless baryon density and the dimensionless 
fluid velocity given by ffTB"]) around their equilibrium values: 

p = — = 1 + ap x + a 2 p 2 + a 3 p 3 + . . . (54) 
Po 

v x = — = av xl + a 2 v x2 + <J 3 v x3 + . . . (55) 

C s 

We now insert the expansions (1541) and fl55l) into (152]) and into (153]) . neglect the terms 
proportional to <r n for n > 2 and organize the equations as series in powers of a and a 2 , 
finding the following equations: 

J dpi dv xl \ 2 J <9p 2 9pi 9pi , #Pi \ n 

i~^ + ^n + \"^ + ^ + ^ + ^r +Pi ^/ =0 (56) 

and 

. / SV 2 Po \ dpi 1 2 / dv xl dv xl dv xl 



( 9v 2 Po \ dpi f g v 2 p \ dp 2 1 (~ 4 ~ VVi ] ( z 7) 

\Mc s 2 m v 2 ) Pl di \Mc s 2 m v 2 J d£ Mp c s L V 3 7 d£ 2 J 1 ' 

Each bracket in the last two equations must vanish: <r{. . .} = and a 2 {. . .} = . The 
final step of RPM is collect the equations from each a order to obtain the nonlinear wave 
equation. From the terms proportional to a we have, in a simple way, after the integration 
with respect to £ and setting the integration constant equals to zero, the following results: 



Vxi = Pi (58) 



and the algebraic relation: 



which gives the sound speed (c s ): 



9v 2 po 



Mc s 2 mv 2 



Mm v 2 



(59) 



From the terms proportional to a 2 we find: 



dp 2 dv x2 dp! dpi dv x i 

-W-^ = ^ + Vxl ^ + pl ^f (61) 

and 

9v 2 Po \dp 2 dv x2 dv xl dv xl dv xl 



Vxi-^r- + Pr 



Mc 2 m v 2 ) di di dr 1 <9£ ' <9£ 
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Qv 2 Po \ dpi 1 / - 4 \ d 2 v x i 

V£ + 77 7 (t+oV) "dn (62) 



Mc s 2 m v 2 J ri d£ Mp c s Ly 3'J d£ 
Using (158]) in ( IBTj) we find: 

^~^f = ^ + 2p1 ^ (63) 
and using (158|) and fl59|) in f )62|) we find: 

^ + 77777 ^ + ^ art ( 64 ) 



<9£ <9£ <9r r <9£ Mp c s L v 3 7 <9£ 
Combining f )63|) and (I64p we find the Burgers equation in the (£, r) space: 

tyi 3 <9 Pl 1 / ~ 4 \ d 2 Pl 

57 + 2* flf = 2Mp^L (, C + 3^ 5e (65) 
Returning (IBTdI) to the (x, t) space (with the help of floU]) and ( lo"Tj) ) we obtain the following 
Burgers equation: 

d d 3 d I f A \ d 2 

where ( l54"j) : <5p# = cpi is the first deviation from the background baryon density po- 



C. Nonlinear waves in a QGP 
From ffTTl) we have: 



(67) 



In the non-relativistic limit e + p = p and from (fl~6j) and (fTTj) we find: 



P^f?) 7/3 W (68) 



3 V2/ 

Inserting ( 1671) and ( )68l) into (|3]) we obtain the following Navier-Stokes equation for the plasma 
phase: 



4 /3\ 7 /3 
3 \2/ 



I 7T 2 / 3 



4/af dv x . duA 4/3\ 7 / 3 2/3 i/ 3 dp B f 4 Nd 2 ^ 

Repeating all the steps shown in the last section, we obtain from the terms of order a: 

v xl = pi and c 2 = 1/3 (70) 

and from terms of order a 2 : 

4 f 3 V /3 ,2/3 n 4/3 <9p 2 4/3\ 7 / 3 2 /3 4/3^x2 16 / 3 \ 7 / 3 2/3 4/3 
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27c, 2 \2/ 



3\ 7 /3 



2/3 „ 4/3 „ 

7r 1 pa 1 pi 



w 1 



^2/3^ 4/3 
7T f p ' 



dv xl 4/3x7/3 
<9r 3 V2/ 



_2/3^ 4/3,, ^1 



dt- 



+ 



? 4 \ d 2 v xl 



(71) 



Inserting (17U|) into (ITTI) and combining this result with ( 163]) we find the following Burgers 
equation in the (x, t) space: 







5pi 



9 A 
c s —dpB 



c s 5p B tt$Pb 
ox 



|(f 



<9t <9x 
where, again, 5ps = op\. 

Comparing the above expression with ( 166]) we see that, apart from small numerical differ- 
ences in the coefficients, nonlinear waves obey the same Burgers equation in both QGP and 
HG phases. We can thus anticipate that, as it was the case of linear waves, the most impor- 
tant differences of nonlinear wave propagation will come from the very different numerical 
values of the viscosity coefficients r\ and (. 



D. Shock waves 



The Burgers equations ( 166]) and ( 172]) may be written in the general form: 



Of 



5pi 



d_ 

' dx 







d 2 



5p B + a 8p B tt^pb = p Tr^fip 



dx 



dx' 



B 



(73) 



where a and p, are the respective nonlinear and dissipative coefficients for the the plasma 
and hadron phase given respectively by: 



a 



OtQGP = C s 
&HG = §C S 



and 



/' 



I^HG ~ A HG 

7/3 



Pqgp 



3 V 2 / 



7T 2 / 3 p 4/3 



x A 



(74) 



QGP 



where for simplicity we have defined the "effective viscosity coefficient" 



A 



AhG = { ChG + \f]HG 



Aqgp = ( (qgp + iVQGP 



(75) 



which assumes distinct values for QGP and HG phases. 



Applying the hyperbolic tangent function method as described in 4lM43| we obtain the 
exact traveling wave solution of (1731) with the free parameter A: 



Sp B {x,t;p:) 



2fi\ 2fi\ 



tanhi X 



o 



a 



a 



x - I c s + - ) t 



(76) 



15 



For the particular choice A = —a/6fi , to ensure 5ps < 1 , the following wave: 



. . . 1 1 , f a 
dp B (x, t;n) = - + -tanh< - — 
3 3 [ 6/i 



x 



a , 
c s + - I 



with the supersonic speed v s 



We follow the analysis performed in 



a 

c s + — 



(77) 



(78) 
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43| of the hyperbolic tangent profile f!77|) . We 



fix t for the QGP and HG phases and obtain the wave profiles shown in Fig. [T] with the 
parameters of Table [H 



TABLE I: Parameters of QGP and HG for Fig. [TJ 





Po (fm 3 ) 


Cs 2 


A (GeVy/m 2 ) 


t(fm) 


QGP 


Wpn 


1/3 


0.05 


0.05 


HG 


1.5 p n 


1/3 


2 


0.05 




FIG. 1: Burgers shock wave solution with different values of the effective viscosity coefficient A for 
QGP and HG. 

In Fig. [JJ we observe that in QGP phase the limit of small /i, in ( 177)) is given by the 
classical shock wave: 

{2/3 for x/t < c s + a/3 
(79) 
for x/t > c s + a/3 
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4l| we rewrite ( 173]) as the following inhomogeneous continuity equation: 

d d ( . a r ,\ d 2 . . . 

di 5pB + Tx [ Cs 5pB + 2 5pB ) =p d^ 5pB (80) 



As in 



where the limiting process p — > given by ( 1791) is a dissipative regularization of the conser- 
vation law: 

_ (5pB + _^^ + _ W j =0 (81) 

For the HG phase in Fig. [I] the baryon perturbation (177)1 is quite different from QGP phase 
and we have: 

5p B (x, t ; pug » Pqgp) = ^ (82) 

due to higher background baryon density and higher viscous coefficients. From these figures 
we can conclude that, due to the large viscosity, shock waves can not be formed in HG. In 
contrast, they can be easily formed in a QGP. This difference may have phenomenological 
consequences. 



E. Wave packets 

The most important feature of the nonlinear equation is that there is a competition be- 
tween the nonlinear and dissipative coefficients. For some particular parameter choices there 
will be a balance between the two terms, in which case localized soliton-like configurations 
may propagate for long distances. The numbers employed here suggest that these localized 
configurations can only survive in the plasma. In order to illustrate this we will perform the 
numerical study of the time evolution of a soliton-like pulse. In this case the initial condition 
is given by: 

SpB(x,t ) = A sech 2 — (83) 
-B . 

for both fl66l) and (1721) . with the amplitude A and width B as free parameters. The viscosity 
coefficients are extrapolated from p, |7[ and represented in the effective viscosity coefficient 
f J75|) . The speed of sound in the HG is calculated by fl60|) and is also c 2 = 1/3 . 

In the analysis of the numerical solutions we must carry out a comprehensive study 
varying each parameter independently, keeping all others fixed. In particular, we study how 



the solution changes when we change the viscosity. Such a study was performed in [23 1 
for the numerical solutions of a similar wave equation, derived from the relativistic Navier- 
Stokes equation applied to the QGP at high temperatures. The conclusion was that the 
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TABLE II: Parameters of QGP and HG for Fig. 





Po (fm 3 ) 


Cs 2 


A {GeV/fm 2 ) 


B(fm) 


QGP 


Wpn 


1/3 


0.05 


0.5 


HG 


1.5 p n 


1/3 


2 


0.5 




FIG. 2: Viscosity effects in hadron and plasma phase, 
most relevant parameter was the viscosity. As for the others, we find that they have here 



the same relative importance as in [23| . In Fig. [2] we compare the hadron gas (on the 
two right panels) with the quark gluon plasma (on the two left panels) for the same scales. 
Comparing left and right we can observe the strong damping which happens in the hadron 
gas, which comes ultimately from the larger values of the viscosity coefficients. Comparing 
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the upper with the lower panels we see that increasing the initial amplitude of the pulse in 
both phases does not change much its evolution, damping and slowing down. In the plasma 



phase we can also observe the formation of a "wall" as can be seen in Fig. 2(a) and Fig. 
[2R| 



V. CONCLUSIONS 



In this work we have studied wave propagation in cold and dense matter, both in a hadron 
gas and in a quark gluon plasma phase. In deriving wave equations from the equations of 
hydrodynamics, we have considered both smaller and larger amplitude waves. The former 
were treated with the linearization approximation while the latter were treated with the 
reductive perturbation method. 

Linear waves were obtained by solving an inhomogeneous viscous wave equation and they 
have the familiar form of sinusoidal travelling waves multiplyed by an exponential damping 
factor, which depends on the viscosity coefficients. Since these coefficients differ by two 
orders of magnitude, even without any numerical calculation we can conclude that, apart 
from extremely special parameter choices, in contrast to the quark gluon plasma there will 
be no linear wave propagation in a hadron gas. 

Nonlinear waves were obtained by solving the Burgers equation, which was derived from 
the equations of hydrodynamics with the RPM. We could find an analytical solution of 
the Burgers equation which describes a shock wave. Varying the viscosity parameters in the 
appropriate range we could conclude that, with the exception of extremely special parameter 
choices, there is no shock formation in a hadron gas. If some external agent would try to 
create a sharp density discontinuity in a hadronic medium, viscosity would immediately 
smooth it and wash it out. In a quark gluon plasma, on the other hand, shocks may be 
formed. Exactly the same features were observed in the wave packet numerical solutions 
of the Burgers equation. In the hadron gas, the wave packet moves very slowly and its 
amplitude is very quickly reduced. In the quark gluon plasma there is a balance between 
nonlinearity and dissipation that prevents the wave from breaking and from dispersion. 
Density perturbations may "survive" longer in QGP than in HG. 

We believe that our work may help in discriminating the quark gluon plasma from the 
hadron gas. The next and most difficult step is to connect the waves studied here with 
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observables and plug our analytical studies in realistic numerical simulations of low energy 
heavy ion collisions or compact stars. Work in this direction is in progress. 
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